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f--^ This paper presents an analysis of n~Pb — » X~Pb — » n~ n~ zr+Pb events at 190 GeV/ c beam 

momentum and very low four-momentum transfer t' < 0.001 GeV 2 /c 2 . Coherent scattering 
. off the nucleus as a whole dominates with contributions from Reggeon, Pomeron and photon 

exchange. The latter originates from Primakoff reactions and is identified by the sharp Coulomb 
peak of intensities at t' « 0. The partial-wave analysis of these data focusses on new techniques 
for the extraction of the Primakoff contribution at low masses. Its measured absolute cross- 
section at \/s <5m n is well in agreement with the prediction from chiral perturbation theory. 
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Oh 1 Meson Spectroscopy at Low Momentum Transfer 

Dissociation of pions on nuclear or hydrogen targets provides clean access to the light 
t— H meson spectrum. In case of heavy target nuclei and low momentum transfer, mesons can 

be produced by two interaction mechanisms: diffractive production by Reggeon f'-channel 
exchange and Primakoff (or Coulomb) production via the exchange of quasi-real photons. 
The latter has its dominant contribution at vanishing momentum transfer and introduces 
spin-projection M = ±1 to the produced system. At low masses, z'.e.at the threshold of 
the investigated final state, the only involved hadrons are pions scattering with each other 
t-h and the photon. Thus predictions from chiral perturbation theory (ChPT), which take into 

account coupling to a real photon simultaneously, are supposed to be applicable. At higher 
masses, resonances can be produced, so that the radiative coupling of the #2(1320), and 
^ potentially also of heavier mesons, can be investigated. Also the interference between the 

Coulomb and diffractive production of these resonances can be studied. 

COMPASS is a multi-purpose fixed-target experiment at the CERN SPS, that investigates 
the structure and spectroscopy of hadrons. The two-stage high-precision spectrometer [1] 
can detect outgoing particles within a large range of scattering angles and particle momenta, 
and provides a uniform acceptance, especially for reactions featuring low to intermediate 
momentum transfer from the beam to the target. For a short run with a 190 GeV/c hadron 
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beam (composed of 96.8% , 2.4% K" and 0.8 %p) on thin lead disk targets in 2004, 
additional silicon micro-strip detectors were installed downstream of the target to resolve 
also smallest scattering angles. With this setup, about 4 million exclusive n~ n~ n + events 
have been collected with a dedicated multiplicity trigger for charged-particle final states. 
About 1 million of these feature very low momentum transfer t' < 0.001 GeV 2 / c 2 , with 
t' = \t\ - \t\ min constituted of the squared four-momentum transfer t from the beam to the 
produced system, and |f |min the minimum value of |f| allowed by kinematics. In this region, 
referred to as "Primakoff t"', photon exchange appears to be relevant in the data. 

Partial-wave analysis (PWA) strives for determining all resonances present in a given data 
set and their properties by fitting angular distributions and taking into account interfer- 
ence effects. The PWA presented here is based on the commonly-known isobar model, 
that was adapted to the specific needs of the analysis. According to this model, the pro- 
duced resonances decay via intermediate two-particle decays into the particles observed 
in the spectrometer. For 3n events, a partial wave in the reflectivity basis is written as 
J PC M e [isobar n]L, defining the quantum numbers of the resonance, } PC , spin projection M, 
reflectivity e, the isobar, and the angular momentum L between the di-pion resonance and 
the unpaired pion. The PWA method and its basic assumptions are described in more detail 
in Ref. [2] and the references therein. 



2 Diagrams from ChPT 



In the Primakoff t' region, decay amplitudes with M = 1 can be attributed to Primakoff 
production, z'.e.effectively 7r~7 — > X~ — > n~n~n + scattering when observing 7r~Pb — > 
X Pb — > n~ n~ n + ¥h. In addition the low mass region, i.e. m^ n < 0.72 GeV/c 2 , is governed 
by pion-scattering instead of resonance production so that ChPT calculations are applicable. 
The leading order processes in ChPT are calculated in [3] and depicted in Fig. 1, including 
also the lead nucleus which is not part of the calculation. As it has been shown in [4] that 
for low masses the next order in ChPT (i.e. loops) does not change the cross-section of wy- 




Figure 1 : Leading order processes in ChPT [3] for the reaction n 7 — > n n n + , embed- 
ded in the Primakoff reaction contribution to 7T~Pb — > n~ n~ 7T + Pb. 
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Figure 2: t' dependences of the M = total intensity (left) and the intensity of the chiral 
amplitude (right) confirm the diffractive production of the M = states as well as the 
Primakoff nature of the intensity picked up by the chiral amplitude. 



scattering, <7 7 , appreciably, it is sufficient to introduce the leading order calculations for a 7 
into the PWA, multiplied with the quasi-real photon density provided by a nucleus of charge 
Z as given by Weizsacker and Williams [5]. The experimentally observed cross-section is 
given by 

da Fh = ci-Z 2 t' ■ F 2 (Q da 7 
dsdt' d<$ n 7r(s - ml) ' (f + t min ) 2 ' d<$>„ 

with d<3>„ the phase-space element for the final state system X~ and F 2 f (t') the effective 
lead form factor, and is introduced in its fully differential form as an additional decay 
amplitude with M = 1 to the PWA fit. This special partial wave, called "chiral amplitude", 
is mathematically not orthogonal to the M = 1 waves of the isobar model, so that these 
have in fact to be replaced by the chiral amplitude in the low mass region. The resulting 
fit quality is compatible with the quality of a previous fit that used up to 6 waves with the 
(7T7r) s wave or the p(770) as isobars, but the chiral-amplitude fit works with much less 
fitting parameters. At masses above m-$ n > 0.72 GeV/ c 2 isobaric decays that compete with 
the chiral amplitude have to be taken into account again. 

In order to get a confirmation of the production mechanism that is attributed to the chiral 
amplitude, a dedicated set of PWA fits was carried out in several bins of t' in the low- 
mass region without taking into account any t' dependence of the decay amplitudes. The 
resulting intensity spectra of the total intensity of the chiral amplitude and the total M = 
intensity in this mass region, dependent on t' , are shown in Fig. 2. Both spectra are fitted by 
an exponential <x exp(— bt') for t' < 0.001 GeV 2 /c 2 . The M = intensity features a slope 
b ~ 400GeV~ 2 /c~ 2 which is compatible with diffractive production off lead nuclei. The 
ChPT intensity is much steeper with b w 1600GeV~ 2 /c~ 2 which is consistent with the 
expected effective slope of Primakoff production if the t' dependence from Eq. 1 is smeared 
with the experimental resolution. 
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Figure 3: Luminosity determination from the free decay of beam kaons: Reconstructed 
invariant mass spectra (left) and decay vertex positions (right) for both COMPASS data 
(top row) and simulated kaon decays in the COMPASS target region (bottom row). 



3 Absolute Cross-Section and ChPT Prediction 

For the comparison with the ChPT prediction from [3], PWA fits for 0.44 < < 
0.72GeV/c were performed in 40MeV/c mass bins for t' < 0.001 GeV 2 /c 2 , that use 
only the chiral amplitude associated with Primakoff production as M = 1 contribution. 

The total intensity of the chiral amplitude, fitted from the experimental data and corrected 
for acceptance effects, is converted into its absolute cross-section via the well-known 
thickness of the lead target and the incoming beam flux. The absolute beam intensity was 
not monitored reliably, but the effective beam flux can be obtained with good precision 
from the observed K~ — > n~ n~ n + decays in the COMPASS target region, which originate 
from the kaon component present in the negative hadron beam. Thus the normalization 
is available from the same data- taking set (z'.e.in the same final state, but with broader 
cut on the decay vertex position) featuring the same systematics concerning the relevant 
trigger and detector efficiencies. Fig. 3 presents the invariant 3zr mass spectrum, obtained 
with broad vertex cut, showing the clean kaon signal (top left), from which the observed 
number of kaons is obtained after subtraction of the small background that is estimated 
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Figure 4: Measured cross-section of n Pb 
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by a linear fit. The corresponding decay vertex distribution (Fig. 3 top right) shows that 
these decay vertices are reconstructed in free space, at the downstream end of the spectrum 
sharply limited by the position of the respective charged-particle multiplicity trigger, and 
upstream limited by the positions of the beam telescope detectors defining the vertex 
reconstruction efficiency. The background from pion interactions in the lead target is 
obtained from the vertex distribution in the neighboring mass region, scaled accordingly 
and subtracted statistically. The selected region of reconstructed vertex positions, that 
proceeds to the mass spectrum discussed before, is selected as indicated to assure uniform 
reconstruction efficiency inside the considered region. These two distributions were also 
confirmed by a dedicated full Monte Carlo (MC) simulation of kaon decays in the relevant 
region of the COMPASS spectrometer. Fig. 3 (bottom left) shows the mass distribution 
of the reconstructed MC kaon decays. The contribution from pion interactions present 
in the real data as smoothly rising background is now completely absent. The shape of 
the reconstructed kaon mass is precisely modeled. The broad part at the base of the kaon 
peak is traced back to kaons decaying upstream of the lead target. In that case all three 
decay pions undergo multiple scattering in the lead, which leads to the observed distortion. 
Fig. 3 (bottom right) presents the distributions for both simulated and reconstructed kaon 
decay vertices, confirming the indicated cuts on the reconstructed vertices in real data to be 
reasonably chosen. 

The black points in Fig. 4 (left) show the cross-section of the Primakoff contribution to 
7T~Pb — > tC tt~ 7T + Pb, i.e.the total intensity of the chiral amplitude in the respective mass 
bins converted by the kaon normalization, with their statistical uncertainties from the PWA. 
The triangles picture the ChPT prediction as given by Eq. 1 which is integrated in the same 
mass bins as used in the analysis of the data. The full systematic uncertainties include the 
uncertainty of the luminosity determination from the kaon decays, results from different 
fitting models, and estimations for the radiative corrections to be applied [6], so that the 
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total uncertainties sum up to about 20%. 

From this measured cross-section the quasi-real photon density is divided out to get the 
cross-section of 7i~y —> n~n~n + presented in Fig. 4 (right). The ChPT calculation from [3] 
is depicted here as a continuous line, while the data points are still placed at the center of 
the respective bin in which they have been obtained, showing good agreement with the 
theoretical prediction. 

4 Conclusion and Outlook 

The result presented in this paper confirms the ChPT leading order calculation for 7r~7 —> 
n~ n~ n + with an experimental uncertainty of about 20%, both in its integrated strength 
and its shape in the 5 dimensional three-body phase space for m^, n < 0.72 GeV/c 2 , for the 
first time. 

For the adjacent region of higher masses, both loop effects and explicit p contributions have 
additionally to be taken into account for the chiral amplitude. The need to include isobaric 
decays with M = 1 becomes obvious for masses approaching the 02(1320) resonances at 
latest. These extensions of the presented analysis are subject of work currently being carried 
out. 
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